function dstate = TerminalDynamic(state, alpha, bank)
global mass S g
z = state(3);
V = state(4);
fpa = state(5);
azi = state(6);
[L, D] = getLandD(z, V, alpha, S);

xdot = V.*cos(fpa).*sin(azi);
ydot = V.*cos(fpa).*cos(azi);
zdot = V.*sin(fpa);
speeddot = -D./mass - g.*sin(fpa);
fpadot   = L.*cos(bank)./mass./V - g.*cos(fpa)./V;
azidot   = L.*sin(bank)./cos(fpa)./V./mass;
dstate = [xdot  ydot  zdot  speeddot  fpadot  azidot];
end